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High-precision mass measurements of more than thirty neutron-rich nuclides around the 
Z=28 closed proton shell were performed with the triple-trap mass spectrometer ISOLTRAP at 
ISOLDE/CERN to address the question of a possible neutron shell closure at A'^=40. The results, 

for 57.60.64-69j^i ^ 38), ^^-^^'^^Cu {Z = 29), and 63-65,68-78q^ ^ 3^)^ j^^^^ ^ ^^,^^5^^ 

certainty of the order of 10~*. In particular, the masses of ^^"^'''^^Cu have been measured for the 
first time. We analyse the resulting mass surface for signs of magicity, comparing the behavior of 
JV=40 to that of known magic numbers and to mid-shell behavior. Contrary to nuclear spectroscopy 
studies, no indications of a shell or sub-shell closure are found for N=40. 



PACS numbers: 21.10.Dr, 21.60.Cs, 27.50. +e, 32.10.Bi 

I. INTRODUCTION 

A striking parallel between the atomic and nuclear sys- 
tems is the occurrence of closed shells. The behavior 
of the atomic system is largely governed by what can 
be considered as an infinitely massive and point-like nu- 
cleus. Describing nuclear behavior, however, is a par- 
ticularly difHcult task given its composition of neutrons 
and protons, similar in mass yet different in charge. The 
nucleon interaction is so complicated that ground-state 
properties are not globally predicted with particularly 
good precision. A property crucial to the understanding 
of the nuclear system is the behavior of its shell structure 
as a function of the varying composition of protons and 
neutrons. The fact that shell structure seems to be mod- 
ified in systems where the number of neutrons TV and the 
number of protons Z are unbalanced {i.e. far from the 
equilibrium region of stable nuclides) is one of the key 
questions of today's nuclear physics research. 

Over the last 20 years, magic numbers have been found 
to vanish in certain region of the chart of nuclides, the 
first one being TV = 20 for sodium H and later, magne- 
sium More recently, = 8 and = 28 0, H] 
have also disappeared. Conversely, "new" magic numbers 
such as A' = 16 [a] and A = 32 0, i, [§| have also been 
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found. One case of particular interest is that of A^ = 40 
because of the unexpected events that have transpired 
since the first studies in 1982. At that time, Bernas 
et al. fiol] showed that the first excited state of 28^140 
was 0+, establishing a new case of 2+ and 0+ inversion. 
This was compared to the case of |QCa2o, a doubly-magic 
nucHde [llj where such an inversion was known. Conse- 
quently, Bernas et al. concluded ^^Ni to be doubly-magic. 

In 1995, Broda et al. [12] published a comprehensive 
summary of spectroscopy work since 1982 and elaborated 
the excited spectrum of ^^Ni, finding the first excited 
state to be 0+ (as Bernas et aL[lO]), 2+ as the second 
excited state and a 5~ isomeric state. As this is the 
same situation for the ^°Zr excited states, they concluded 
that ^®Ni was spherical, implying a significant sub-shell 
closure at A^ = 40. Shell-model predictions of isomeric 
states near magic nuclides motivated the experimental 
investigations of Grzywacz et al. (l3l | in 1998. They dis- 
covered many isomeric states in the vicinity of ^^Ni, fur- 
ther strengthening the case for its doubly-magic charac- 
ter. In 1999, /3-decay studies were carried out by Han- 
nawald et al. ^14j] . who found long half-lives for the neigh- 
boring isotones (copper, manganese) at A^ = 40 indicat- 
ing an increase in collectivity. However, /3-decay studies 
by Mueller et al. (l5l | the same year showed that the sta- 
bilizing effect of A^ = 40 disappeared when moving away 
from 68Ni. 

The powerful tool of Coulomb excitation was brought 
to bear on ^*Ni in 2002 when SorHn et al. |16|] measured 
the B[E2) value (which is the probability of transition 
between the ground state O"*" and the excited state 2+). 
B{E2) is expected to be small for magic nucHdes which 
are difficult to excite, and to be large for deformed nu- 
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elides. The measured B{E2) value was unexpectedly 
small, reinforcing the magic nature of ^^Ni. Sorlin et al. 
attributed the lack of corroborating evidence from the 
mass surface to an erosion of the TV = 40 sub-shell, ero- 
sion confirmed by recent measurements [13, How- 
ever, a concerted theoretical effort published by Lan- 
ganke et al. [31 argued against the doubly-magic nature 
of ^^Ni, noting that the "missing" B{E2) strength lies at 
much higher energy (>4MeV). 

According to Bohr and Mottelson [IS]: "In terms of the 
expansion of the total binding energy, the shell structure 
appears as a small correction compared to the surface 
energy... Despite the smallness of these effects on the 
scale of the total nuclear energy, they are of decisive 
importance for the structure of the low-energy nuclear 
spectra..." In the Hght of these conflicting experimental 
and theoretical signatures as well as the relatively large 
uncertainty on the binding energies in this interesting re- 
gion, high-precision mass measurements were carried out 
with the mass spectrometer ISOLTRAP in an attempt to 
bring some clariflcation to this situation. Time-of-flight 
mass measurements had been performed in 1994 [2lj but 
although they gave no indication that TV = 40 was magic, 
the precision was insufficient to be conclusive. The most 
accurate mass measurements today are performed in Pen- 
ning traps [H and ISOLTRAP at CERN has pio- 
neered the application to radioactive nuclides [zi. |25| . 
The experimental setup of ISOLTRAP is presented in 
sectionlTTl and the measurements in the region of iV = 40 
and their evaluation are described in section lllll A com- 
parison to mass models follows in section IIVI and the 
question of = 40 is discussed in the light of the new 
results in the last section. 



II. THE ISOLTRAP SETUP 

A. Experimental setup 

ISOLTRAP is a high-precision Penning-trap mass 
spectrometer, located at CERN's ISOLDE facility [H 
which delivers mass-separated beams of radionuclides. 
ISOLTRAP is composed of three main parts (see Fig.[T]). 
First, a linear gas-filled radio-frequency quadrupole 
(RFQ) trap, used as cooler and buncher, adapts the 60- 
keV ISOLDE ion beam to the ISOLTRAP requirements 
with respect to kinetic energy, time structure, and beam 
emittance [131 ■ The second part is a gas-filled, cyHndrical 
Penning trap [H] in which a mass-selective helium buffer- 
gas cooling technique [22(| with a resolving power of up to 
10^ is used for isobaric cleaning. This preparation trap is 
installed in a 5=4.7 T superconducting magnet. Finally, 
the cooled ion bunch is transferred to the precision Pen- 
ning trap for isomeric separation (when required) and 
mass measurement. The precision Penning trap is in- 
stalled in a second superconducting magnet {B = 5.9T). 
The mass is determined by measuring the true cyclotron 
frequency Vc = qB/{2Trm) of the stored ion (see next 



paragraph). The magnetic field B is determined from a 
measurement of the cyclotron frequency of a reference ion 
whose mass is well known. The setup also includes an off- 
line ion source to produce stable ions, used as reference 
masses. 



B. Mass measurement procedure 

Ion confinement in a Penning trap is based on the ap- 
plication of an electrostatic field and a magnetic field to 
store ions in the axial and radial directions, respectively. 
The ion motion in a Penning trap is a superposition of 
three independent harmonic oscillator modes, one in the 
axial direction with frequency ly'z and two in the radial di- 
rection, i.e. the cyclotron motion with reduced frequency 
and the magnetron motion with frequency ly- f3d,'3l|. 
In a purely quadrupolar electric field, the frequencies are 
related as follows: 

Vc = iy+ + i^-. (1) 

Ion beams are alternatively delivered from ISOLDE or 
from an off-line ion source and injected into the RFQ, 
mounted on a 60-keV pedestal, where they are cooled 
and bunched. The ion bunch from the RFQ is sent to 
the preparation trap. Ion collisions with the buffer gas 
inside this trap first cool the axial motion. A dipolar ex- 
citation with a frequency v- is then appHed to increase 
the magnetron radius of all ion species, making it larger 
than the exit hole of the trap. To select the ions of in- 
terest, an azimuthal quadrupole radio-frequency electric 
field at frequency is appHed which couples the radial 
modes. Since one mode is cooled by the gas, the radius 
is reduced and the ion cloud is centered. In this way the 
trap works as an isobar separator with a resolving power 
R = m/Am of 10'' to 10^^ [H. 

The purified ion beam is transferred to the precision 
trap, where different excitations are performed. A phase- 
sensitive dipolar excitation at is applied to increase 
the magnetron radius of the ion motion [33|. If there 
are contaminants (isobars or isomers), a second, mass- 
dependent dipolar excitation is performed at v+ to re- 
move them |33|- Finally, an azimuthal quadrupole radio- 
frequency field is applied to convert the initial magnetron 
motion into cyclotron motion. At i^fiF = Vc, a full con- 
version is obtained, leading to an increase of the orbital 
magnetic moment /x and the associated radial kinetic en- 
ergy E = /iSfs^l- After ejection at low axial energy, 
ions pass the inhomogeneous part of the magnetic field 
on their way to an MCP detector (recently replaced by a 
channeltron detector [s^) at the top of the setup. Since 
the axial acceleration in this fringe field is proportional to 
/i • dB/dz, the shortest time of fiight (TOF) is observed 
for vrf = i^c [sil. 

The mass resolution in the precision trap depends 
strongly on the conversion time used for the excitation. 
The line width Av of the resonance is mainly determined 
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FIG. 1: Sketch of the experimental setup of the ISOLTRAP mass spectrometer, including the main parts: a gas-filled linear 
radio-frequency quadrupole (RFQ) trap for capturing and preparing the ISOLDE beam, a gas- filled cylindrical Penning trap for 
isoharic separation, and a hyperbolic Penning trap for the mass measurement. The micro-channel plate (MCP) detectors are 
used to monitor the ion transfer and to measure the extracted-ion time of flight (TOP) together with the channeltron detector. 
The inset presents a time-of-flight (TOF) cyclotron resonance for radioactive Ni^ ions. 



by the duration of the appHed RF-field (Trf) used to 
couple the two radial motions. The relation is |34|: 



Aiy{FWHM) 



0.9 

Trf 



(2) 



The statistical precision in the cyclotron frequency de- 
termination is given by [37i |: 



oc 



(3) 



with N being the number of ions and R = vTrf the re- 



solving power. With sufficiently long excitation times 
(few seconds), a resolving power of up to 10'' can be 
reached. As an example of a cyclotron frequency mea- 
surement, the inset of Fig.[l] presents the time-of-flight 
(TOF)-resonance curve of one of the two measurements 
of radioactive ^^Ni. The mean TOF of the ions as a 
function of the applied radio-frequency (RF) is shown. 
The solid line is a flt of the well-known line-shape [31| to 
the data points. This measurement was performed with 
about 1000 ions using an excitation time Trf — 900 ms, 
resulting in a resolving power of 1.1 x 10^ and a relative 
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frequency uncertainty of ^vjv = 6 x 10 



III. MEASUREMENTS OF THE NI, CU, AND 
GA ISOTOPES 

The nuclides 57,60,64-69]N^j^ 65-74,76 (-.^^^ ajjj 

63-65,68-78Qg^ have been investigated with ISOLTRAP. 
They were produced at ISOLDE by bombarding a 
uranium carbide (UC) target with 1.4-GeV protons from 
CERN's Proton Synchroton Booster. The ionization 
was achieved for gallium with a tungsten (W) surface 
ionization ion source and for copper and nickel with 
the resonance ionization laser ion-source (RILIS) [s^ . 
ISOLDE'S General Purpose Separator (GPS), with a 
mass resolving power of about 1000 was used. The 
proton- rich isotopes ^■^~^^Ga were measured in a differ- 
ent experiment using a ZrO target and ISOLDE'S High 
Resolution Separator (HRS), which has a mass-resolving 
power of about 3000. Both targets were bombarded 
using pulses containing up to 3 x 10^^ protons. 

The yields of nickel and copper were fairly intense at 
about 10^ ions/s. The efficiency of ISOLTRAP is better 
than 1% so a beam gate was used in order to limit the 
number of ions sent to the precision trap and minimize 
ion-ion interactions that cause frequency shifts. The typ- 
ical number of ions simultaneously stored in the precision 
trap was between 1 and 8. 

Despite the good yields of nickel and copper nuclides, 
up to three orders of magnitude more surface-ionized gal- 
lium was present. For the measurement of ^*Ni shown in 
Fig.[Tl a cleaning of ^*Ga was applied in the prepara- 
tion trap. The ratio between the yield of ^^Ga and ^^Ni 
was "only" a factor of ten which was low enough to al- 
low an effective cleaning. This ratio was higher farther 
from stability and prevented the measurement of more 
neutron-rich nickel and copper since the preparation trap 
was saturated by the gallium isobars. Similarly, a sig- 
nificant contamination of titanium oxide prevented the 
measurement of more proton- rich gallium isotopes, and 
the presence of rubidium isobars made the measurement 
of more neutron-rich gallium isotopes impossible. 

The results from the data analysis is the ratio 
Vc,ref /i^c [1^, since the atomic mass m of the ions is cal- 
culated from the ratio between the cyclotron frequency 
of the reference ion Vcref and the cyclotron frequency of 
the ion of interest Vc, the atomic mass of the reference 
^^Rb fii^l, and the electron mass 

m= [mgsRb - me) + me. (4) 

All the results were evaluated in order to include them 
in the Atomic-Mass Evalution (AME) table [ilj. The ta- 
ble of atomic masses results from an evaluation of all 
available experimental data on masses, including direct 
measurements as well as decay and reaction studies. The 
AME forms a linked network and uses a least-squares ad- 
justment to derive the atomic masses. Among all nuclear 



ground-state properties, such an evaluation is unique to 
mass measurements. 

The mass values from the present measurements are 
presented in TablesU (Ni), [11] (Cu), and |IlI] (Ga). These 
tables give the ratio of the cyclotron frequency of the 
^^Rb"*" [40?! reference mass to that of the ion of interest. 
The corresponding uncertainty takes into account a sta- 
tistical uncertainty depending on the number of ions, and 
a systematic error [39| . The derived mass excess value 
is indicated for comparison with the AME tables from 
1995 and 2003. Since the latest Atomic-Mass Evaluation 
(AME2003 [Hi) includes the data from this work, the in- 
fiuence of the ISOLTRAP measurements is also provided. 
Among the 36 nucHdes measured here, the influence is 
100% for 22 of them. 

The nickel results are presented in TableU and in 
Fig. [21 This figure presents the difference between the 
mass excess measured by ISOLTRAP and the AME1995 
values. Note that even for the stable nickel isotopes 
the precision of the mass values is improved. With 
the exception of ^^Ni (see below) the results are in 
good agreement with the 1995 table but much more 
precise. The masses of 57,60,65]\jj agree with the 1995 
table within the error bars, and were measured with 
the same order of uncertainty. The combination of 
the previous value and the ISOLTRAP measurement 
reduces the final uncertainty. The results contributing 
to the ^^Ni mass value are presented in Fig.[3l This is a 
special case because it is in strong disagreement with the 
AME1995 table [43]: a difference of more than 400 keV 
was observed. The AME1995 value was derived from 
a ™Zn(i'^C,i^O)S''Ni reaction [Hi and a time-of-fiight 
measurement [2l|. The ISOLTRAP value disagrees 
with the value from the reaction but is in agreement 
with the time-of-fiight measurement. Since the value 
of ISOLTRAP is much more precise, the AME2003 
includes only this value. 



The copper results are Hsted in TablelTTl a compar- 
ison with the AME1995 values is given in Fig.[4l An 
improvement of the mass uncertainty was achieved 
for all investigated copper isotopes. The values are in 
good agreement with previous values, except for ''°Cu". 
This important difference is due to an incorrect state 
assignment. ISOLTRAP's high resolving power of more 
than 10^, in combination with /3-decay studies and 
selective laser ionization allowed us to perform a clear 
identification of each state [45]. Moreover, this high 
resolvin g p ower allowed us to resolve isomeric states 
in ^^Cu [401 and to measure them independently. The 
masses of 72-74,76 

were previously unknown. They 
are compared to model predictions in Section lTVl 
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TABLE L ISOLTRAP results for nickel isotopes: nuclide; half life; frequency ratio Va,re.f/vc of nickel isotope to reference 
nuclide Rb'^ fJ^J, corresponding mass excess (ME); mass excess from AME1995; new mass excess from AME2003; influence 
of the present result on the AME2003 value. 



Isotopes 


Half life 






ISOLTRAP 


AME1995 


AME2003 


Influence on 




2^1/2 






ME (keV) 


ME (keV) 


ME (keV) 


AME2003 


iNl 


oo.D n 


u.o / uo / oooyo 


(31bj 


-00Uo4.z \^.o) 


-ooU/o.o ^.i.yj 


-odUoz.U (i.oj 


oz.U /o 


«°Ni 


Stable 


0.7057986239 


(183) 


-64472.7 (1.4) 


-64468.1 (1.4) 


-64472.1 (0.6) 


16.6% 




Stable 


0.7528734602 


(163) 


-67096.9 (1.3) 


-67095.9 (1.4) 


-67099.3 (0.6) 


21.9% 


65 Ni 


2.5 h 


0.7646753441 


(285) 


-65129.0 (2.3) 


-65122.6 (1.5) 


-65126.1 (0.6) 


7.8% 


66 Ni 


55 h 


0.7764412560 


(181) 


-66006.3 (1.4) 


-66028.7 (16.0) 


-66006.3 (1.4) 


100% 


6^Ni 


21 s 


0.7882468785 


(362) 


-63742.7 (2.9) 


-63742.5 (19.1) 


-63742.7 (2.9) 


100% 


68Ni 


29 s 


0.8000274080 


(377) 


-63463.8 (3.0) 


-63486.0 (16.5) 


-63463.8 (3.0) 


100% 


69 Ni 


12 s 


0.8118484759 


(466) 


-59978.6 (3.7) 


-60380 (140) 


-59979 (4) 


100% 



TABLE II: ISOLTRAP results for copper isotopes: nuclide; half life; frequency ratio i'c,ref/i'c of copper isotope to reference 
nuclide Rb'^ fj^] - corresponding mass excess (ME); mass excess from AME1995; new mass excess from AME2003; influence 
of the present result on the AME2003 value. Previously unknown values derived from systematic trends are marked with #. 



Isotopes" 


Half life 








ISOLTRAP 


AME1995 


AME2003 


Influence 




Ti/2 








ME (keV) 


ME (keV) 


ME (keV) 


on AME2003 




Stable 





7646483448 


(139) 


-67264.5 


(1.1) 


-67259.7 (1.7) 


-67263.7 


(0.7) 


36.8% 




5.1 m 





7764380632 


(257) 


-66258.8 


(2.0) 


-66254.3 (1.7) 


-66258.3 


(0.7) 


11.1% 




62 h 





7882016658 


(155) 


-67318.8 


(1.2) 


-67300.2 (8.1) 


-67318.8 


(1.2) 


100% 




31.1 s 





8000008176 


(199) 


-65567.0 


(1.6) 


-65541.9 (45.6) 


-65567.0 


(1.6) 


100% 




3.7 m 





8000098791 


(188) 


-64850.3 


(1.5) 


-64818 (50) 


-64845.4 


(1.7) 


50% 




2.8 m 





8117756816 


(174) 


-65736.2 


(1.4) 


-65739.9 (8.1) 


-65736.2 


(1.4) 


100% 


™Cu« 


45 s 





8235875816 


(199) 


-62976.1 


(1.6) 


-62960.3 (14.5) 


-62976.1 


(1.6) 


100% 




33 s 





8235888547 


(258) 


-62875.4 


(2.0) 


-62859 (15) 


-62875.4 


(2.0) 


100% 


™Cu" 


6.6 s 





8235906419 


(272) 


-62734.1 


(2.1) 


-62617 (15) 


-62734.1 


(2.1) 


100% 


^^Cu 


19 s 





8353679363 


(194) 


-62711.1 


(1.5) 


-62764.2 (35.2) 


-62711.1 


(1.5) 


100% 


^^Cu 


6.6 s 





8471819597 


(182) 


-59783.0 


(1.4) 


-60060# (200#) 


-59783.0 


(1.4) 


100% 


^^Cu 


4.2 s 





8589690332 


(491) 


-58986.6 


(3.9) 


-59160# (300#) 


-58987 


(4) 


100% 


^^Cu 


1.6 s 





8707837184 


(779) 


-56006.2 


(6.2) 


-55700# (400#) 


-56006 


(6) 


100% 




640 ms 





8944013229 


(843) 


-50976.0 


(6.7) 


-50310# (600#) 


-50976 


(7) 


100% 



"g,m,n denote the ground, first excited, and second excited state, 
respectively, of tlie nuclide. 



The gallium results are presented in Tablelllll and in 
Fig. [SI The ^*Ga mass uncertainty, 6m/ m « 5.4 • 10"'' 
is much higher than for all the other nuclides. This 
is due to the use of a shorter excitation time (100 ms 
compared to 900 ms for the other nuclides) and to a lack 
of statistics: only 530 ions were observed, compared to 
at least 3000 for most of the other ones. The ISOLTRAP 
value is still in agreement with the AME1995 value but 
has no influence. For all other gaUium isotopes measured 
by ISOLTRAP the uncertainty was decreased. For five 
of them, it was decreased by more than a factor of 20, 
and for ^^Ga, almost 100 times. 

The case of ^^Ga was complicated by the possible pres- 



ence of a 9.5-second isomeric state having an excita- 
tion energy of only 60 keV (this accounts for the large 
AME1995 error bar in Fig. [5]). Spectroscopy studies per- 
formed in parallel with the mass measurements revealed 
no indication that the isomer was produced. A two- 
second excitation time was used in order to resolve this 
state in the precision trap but it was not seen. More- 
over, the z-class analysis |33| was performed to examine 
any dependence of the result as a function of ion number, 
but revealed no indication of a contaminant. Therefore 
we are confident that the present result is that of the 
ground-state mass. 
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TABLE III: ISOLTRAP results for gallium isotopes: 
nuclide Rb'^ fjdj, corresponding mass excess (ME); 



nuclide; half life; frequency ratio fc.ref/i'c of gallium isotope to reference 
mass excess from AME1995; new mass excess from AME2003; influence 



Isotopes 


Half life 


In 


IbOLlHAP 


AMi!jl995 


AME2003 


T 

Innuence 




J 1/2 




Mtj (keV) 


Mtj (keV) 


ME (keV) 


on AME200O 


«3Ga 


32 s 


0.7412298391 (167) 


-56547.1 (1.3) 


-56689.3 (100.0) 


-000^ / . i 




100% 


«*Ga 


2.6 m 


0.7529779275 (294) 


-58834.1 (2.3) 


-58834.7 (3.9) 


p;QQQ/j Q 


yZ.K) ) 


75.2% 


^^Ga 


15 m 


0.7647065938 (176) 


-62657.3 (1.4) 


-62652.9 (1.8) 


-OzOO / .z 


yV.o ) 


35.6% 


««Ga 


68 m 


0.799981231 (431) 


-67116.2 (34.1) 


-67082.9 (2.0) 


-0 / UoO. i 


yL.O ) 


0% 


««Ga 


Stable 


0.8117302720 (193) 


-69327.9 (1.5) 


-69320.9 (3.0) 


-oyoz / .o 


I 1 . Z I 


65.3% 


^"Ga 


21 m 


0.8235125549 (272) 


-68910.3 (2.2) 


-68904.7 (3.1) 


-68910.1 


(1.2) 


31.8% 


^^Ga 


Stable 


0.8352740255 (357) 


-70138.9 (2.8) 


-70136.8 (1.8) 


-70140.2 


(1.0) 


13.3% 


^^Ga 


14.1 h 


0.8470706093 (182) 


-68590.2 (1.4) 


-68586.5 (2.0) 


-68589.4 


(1.0) 


53.0% 


^^Ga 


4.8 h 


0.8588335898 (208) 


-69699.4 (1.7) 


-69703.8 (6.3) 


-69699.3 


(1.7) 


100% 


^*Ga 


8.1 m 


0.8706314521 (469) 


-68049.6 (3.7) 


-68054.0 (70.7) 


-68050 


(4) 


100% 


^^Ga 


130 s 
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FIG. 2: Difference between the ISOLTRAP mass excess re- 
sults for nickel isotopes and the AME1995 values lAS] . Dashed 
lines represent the ISOLTRAP error bars. 



FIG. 3: Mass excess of Ni determined by the re- 
action ^° Zn(^'^C}^ Of^ Ni fI2]. and a time-of-fliqht mea- 
surement fs^J, the resulting AME1995 value flsij. and the 
ISOLTRAP value. The AME2003 valueUij differs by 400 
keV with an uncertainty 30 times smaller than the AME1995 
value. 
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FIG. 4: Difference between ISOLTRAP mass-excess values for 
copper isotopes and the 1995 AME values [AS]. Dashed lines 
represent the ISOLTRAP error bars, g denotes ground states 
and m,n isomeric states. 



IV. MASS-MODEL PREDICTIONS COMPARED 
WITH NEW DATA 

Various models and formulae have been developed over 
the years to predict properties of nuclides, particularly 
their mass. A review can be found in [47] where a subset 
of mass models was singled out for comparison. We have 
chosen to compare our experimental data to those, as 
described below. 

The venerable Bethe-Weizsacker mass formula (48l.l49|. 
was based on the liquid drop model and did not include 
shell effects. The nuclear mass m is given by 

m{N, Z)c^ ZrripC^ + NrrinC^ - a^A + a^A^/^ 

+ a,z2A-V3 + a,^„l^_^, (5) 

where nip and m„ are the proton and neutron masses, 
and A the mass number of the nucleus. The parame- 
ters are: the volume term, Os the surface term, Oc 
the Coulomb parameter, and asym the asymmetry pa- 
rameter. Note that the tabulated masses are those of 
the neutral atoms, not of the bare atomic nuclei. While 
inappropriate for mass predictions, it can play an inter- 
esting diagnostic role concerning closed shell effects (see 
section lV^ . 



FIG. 5: Difference between ISOLTRAP mass-excess values for 
gallium isotopes and the 1995 AME values l^]. Dashed lines 
represent the ISOLTRAP error bars. 

For many years, a hybrid approach was adopted for 
predicting masses based on a combination of the macro- 
scopic liquid drop model and microscopic {e.g. shell) 
corrections. The most developed form of these so- 
called mic-mac models is the Finite Range Droplet Model 
(FRDM) [131. 

The Duflo-Zuker (DZ) mass formula [5l|, is a global 
approach, derived from a Shell-Model Hamiltonian and 
gives the best fit to the known masses. Shell-Model cal- 
culations, while well-suited for excitation energies, are 
less so for mass predictions although some efforts were 
made in this direction jHI. 

In the last few years, Hartree-Fock Bogolioubov (HFB) 
calculations have been appHed to the construction of 
complete mass tables. Skyrme forces have tradition- 
ally aimed at predicting a wide range of nuclear prop- 
erties [H, [13, [H, IH| ■ The first microscopic Skyrme-force 
mass formula HFBCS-1 [13, [H| was rapidly superceded 
by HFB-1 [1^ which, in turn, was considerably revised, 
resulting in HFB-2 [60|. A systematic study of the differ- 
ent adjustable parameters followed, resulting in a series 
of formulas up to HFB-9 [Hi, 62, 63, 6||. 

In addition to DZ and FRDM, the ISOLTRAP results 
are therefore compared to HFB-2 and the recent HFB-8 
(HFB-9 did not change the mass predictions apprecia- 
bly). 

One characterization of a model is the root-mean- 
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square (rms) deviation from the mass values to which 
its parameters were fitted, defined by 



N 



\ 



N 



exp 



i \2 
th) ' 



(6) 



where N is the number of experimental rriexp and theo- 
retical rrith masses being compared. A more complete 
description of the rms deviation, including errors, can be 
found infill. Table ITVl shows arms for the models com- 
pared with the AME95 table [i^, which does not include 
the present ISOLTRAP results, and with AMEO3 0, 
which does. Our results improved the overall agreement 
for the HFB models, worsened it for the Duflo-Zuker 
(DZ) mass formula and for FDRM there is no change. 
Examining the isotopic chains individually, we see that 
in all cases the HFB models improved and the DZ 
model worsened. For the FRDM, the better fit for the 
gallium isotopes counters the worse fit for copper and 
nickel. The differences are admittedly small (between 
1 and 10%). While it is tempting to conclude that the 
comparison of the arms might be a demonstration of the 
positive evolution of HFB-2 to HFB-8, it is important to 
recall that unlike FRDM and DZ, HFB-8 was adjusted 
to the masses of the AME03. 




34 36 38 40 42 44 46 48 50 

Neutron number N 



FIG. 6 

model 



Mass 



'erence between ISOLTRAP results and 
ions for the copper isotopes. Note that 
72,73,74,76 measured for the first time and that the more 

recent parameter fit for HFB-8 included these results. 



TABLE IV: The root-mean-square deviation arms (in MeV) 
for different models: the Duflo-Zuker (DZ) mass formula, 
the Finite Range Droplet Model (FRDM), and the Hartree- 
Fock Bogolioubov (HFB) calculations, performed with the 
AME tables of 1995 and 2003 (the latter includes the present 
ISOLTRAP data). Calculations were made for the nickel, 
copper, and gallium isotopes measured by ISOLTRAP. The 
first two rows present the calculation for all nuclides and the 
following rows describe the results for each isotopic chain sep- 
arately. 



Nuclide 


AME Table 


DZ 


FRDM 


HFB-2 


HFB-8 




Ni,Cu,Ga 


AME95 


0.434 


0.555 


0.843 


0.550 


Ni,Cu,Ga 


AME03 


0.451 


0.555 


0.801 


0.530 


Ni 


AME95 


0.623 


0.445 


1.211 


0.732 


Ni 


AME03 


0.640 


0.476 


1.174 


0.678 


Cu 


AME95 


0.426 


0.471 


0.644 


0.601 


Cu 


AME03 


0.451 


0.530 


0.626 


0.563 


Ga 


AME95 


0.280 


0.644 


0.654 


0.375 


Ga 


AME03 


0.291 


0.614 


0.648 


0.384 



Of particular interest for mass models is to compare 
predictions as far as possible from what is already known. 
In the case of the copper isotopes presented here, four 
new masses were determined and one of them (^®Cu) has 



five neutrons more than the most neutron-rich previously 
known mass. The differences of the new ISOLTRAP cop- 
per masses with respect to the above-mentioned models 
are shown in Fig.[6l 

Despite going significantly farther from stability, it is 
difficult to asses which model does a better job. The one 
closest to the new mass of ''^Cu is HFB-8, however the 
other models are not far away. The rras errors on just 
the four previously unknown masses are also similar with 
DZ (0.309 MeV) seeming to follow with a better trend 
compared to all the others (HFB-8: 0.400 MeV; HFB- 
2: 0.566 MeV; FRDM: 0.603 MeV). It is surprising that 
despite all models having their parameters adjusted to 
the mass tables that included those nuclides with < 43, 
those masses are not very well reproduced locally. 

Some nucleon-nucleon effective interactions - like for 
instance Skyrme SKM*, SLy4, or Gogny Dl - are de- 
signed to give rise to a realistic mean field (including pair- 
ing) . They are therefore parameterized on the ground of 
a few available nuclear data for which mean field (includ- 
ing pairing) effects can be reasonably disentangle from 
long range correlations ones (for instance, binding en- 
ergies of doubly magic nuclei only). Such approaches 
of nuclei in which long range correlations are not intro- 
duced in the mean field in an effective and somewhat 
uncontrolled manner do not have as objective to give a 
precise mass formula at the mean (HFB) (including pair- 
ing) level, but to constitute the mean field input of more 
elaborated descriptions of nuclei considering - at least 
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FIG. 7: Difference of the nickel results from the Atomic 
Mass Evaluation 2003 (AME2003) which already includes the 
present ISOLTRAP data and those predicted by HFB-DIS 
(Gogny) and GCM-GOA as a function of neutron number N 
for (left) the mass and (right) the two-neutron separation en- 
ergy. 



some - long range correlations up to the best and there- 
fore able to describe "beyond" mean field a large class 
of nuclear observable (mass formula but also low energy 
spectroscopy, shape coexistence, and transitions, etc ...). 
In this frame, we have performed triaxial HFB calcula- 
tions, using numerical methods and codes described in 
[13, with the Gogny DIS force (sl, [13, IHl . Fig.[7l(left) 
presents the differences between the measured Ni masses 
and those predicted by HFB-DIS, as a function of N. 
There is a large offset {rms difference of 2.473 MeV) for 
the HFB-DIS masses, expected, as explained above, spe- 
cially for mid-shell nuclei where long range correlations 
play an important role. Under these assumptions, we 
could expect at least that the derivative of these quanti- 
ties might be closer to reality. Therefore, in Fig.[7|(right), 
we have plotted the two-neutron separation energy S2n 
[see eq. ((Tj)] derived from the same results. The result is 
encouraging, with an rms deviation of only 0.508 MeV. 

In general, due to the existence of long range correla- 
tions beyond mean field, a unique HFB wave function is 
not well suited to describe the nuclear system. Thus, a 
configuration mixing approach already described and ap- 



plied with some noticeable successes to different nuclear 
problems, for instance to shape coexistence and tran- 
sitions in light mercury isotopes [s^, or Normal-Super- 
deformed phenomena (toI. ItH has been considered. Using 
a Generator Coordinate approach under Gaussian Over- 
lap Approximation (GCM-GOA) in a space constituted 
by HFB (DIS) states under axial and triaxial quadrupole 
constraints allows in this model to treat on the same foot- 
ing rotation and quadrupole vibrations. This approach 
which takes explicitly into account these important cor- 
relations, has been applied to the calculation of nickel 
masses, and the results are shown in Fig. [7] for compari- 
son. Already the mass values (left) are greatly improved 
{rms difference of 0.701 MeV), as are the mass deriva- 
tives (right, rms difference of 0.335 MeV). It would ap- 
pear that going beyond the mean field is to be encouraged 
for future mass predictions. Works in this spirit are also 
underway on the ground of Skyrme forces (see e.g. [t^I). 

V. ANALYSIS OF THE MASS SURFACE 
AROUND Z=29 AND N=40 

As recalled in the introduction, Bohr and Mottel- 
son jlOl explain that the effects of binding energy on nu- 
clear structure are subtle but decisive. As such, accu- 
rate mass measurements are important in order to finely 
analyse the mass surface, notably its derivatives. In this 
section we examine several mass-surface derivatives and 
variations. 



A. Study of the two-neutron separation energy 

The two-neutron separation energy (iS'2n) given by 

S2n{N,Z)=B{N,Z)-B{N -2,Z), (7) 

with B for the binding energy, is remarkable for its reg- 
ularity between shell closures. Generally, S2n decreases 
smoothly with N and shell effects appear as discontinu- 
ities. In the past, discontinuities of S2n versus N were 
often traced to inaccurate Qp endpoint measurements 
and measurements with more reliable, direct techniques 
restored the regularity (see, for example, [t^] for the area 
around ^°®Pb). Hence, part of the motivation was to 
confirm any mass surface irregularities in the TV = 40 
region. Fig. [8] presents the S2n values, from iV = 36 to 
50, prior and after the ISOLTRAP mass measurements. 
Most of the irregularities e.g. at iV = 41 for gallium 
are confirmed. Moreover, the plot reveals a deviation 
from the linear trend between iV = 39 and iV = 41 for 
nickel, copper, and gallium. Also irregularities for gal- 
lium (A^ = 46 — 49) and copper (A^ = 43 — 46) are visible. 

To study the structure more closely we subtract a lin- 
ear function of N determined by the S2n slope preceding 
the purported shell closure. The resulting reduced 5*2,1 
values are presented in Fig. [9] in the region of TV = 82 
(for comparison) and A^ = 40. The A^ = 82 shell closure 
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FIG. 8: Two-neutron separation energies (S-zn) for iron (Z = 26) to germanium (Z = 32) around N = 40. Dashed lines corre- 
spond to the data before the ISOLTRAP measurements. Points with large error bars were not directly measured by ISOLTRAP 
but their value was changed by the link to the measured masses. 



is clearly visible on this plot: there is a change of slope 
between = 82 and N = 84. From these observations 
we can analyse the behavior in the TV = 40 region: there 
is a similar effect between = 39 and iV = 41 where 
the break can be seen at iV = 39 and not at TV — 40, 
surprising for an odd number. The magnitude of this de- 
crease is far smaller (between 500 keV and 1 MeV) than 
the one for the major shell closure at = 82 (around 
4 MeV) . A similar structure is seen between TV = 39 and 
= 41 for nickel, copper, and gallium, but this is not 
an indication of shell closure. It is strange that the same 
structure is visible for both nickel (even Z) and gallium 
(odd Z) whereas germanium is smooth and little is seen 
in the case of zinc. Further measurements to reduce the 
uncertainty on the neighboring cobalt isotopes will be 
needed. 



B. The shell gap 

The neutron shell gap, defined as 

Aiv(TV,Z) = S2niN,Z)~S2n{N + 2,Z) (8) 
= 2B{N, Z) - B{N - 2, Z) - B{N + 2, Z), 

is a good indicator of shell strength. The shell gap def- 
inition is usually only valid for spherical nuclides, i.e. 
around magic numbers. Here, we examine the case of 



A^ = 40 and also investigate how mid-shell gaps compare 
in strength and comportment. Fig.[lOl calculated from 
AME2003 datajUl, shows the shell gap as a function of 
the proton number Z for for various A^. This highlights 
the large shell gap values for magic neutron number with 
peaks at magic Z. It also shows that for A^ = 50 there is 
a peak at Z = 39, and not Z = 40, which is known to be 
semi-magic. This behavior is probably due to the odd- 
even effect in the two-proton separation energy S2p. Not 
surprisingly, the mid-shell-gap (A^ = 39, 66) energies are 
quite small. From this point of view, the case of A^ = 40 
resembles a mid-shell rather than a magic number. 

Fig. E] shows the details of adjacent shell gaps An as 
a function of the proton number Z for different regions: 
(a) around a shell closure, (b) in the region of interest, 
and (c) in a mid-shell region. In Fig.fTlTa). the behavior 
of a strong shell closure is shown for A^ = 82 which is a 
magic number: there is a large difference between N = 82 
and A^ — 81, 83 and the corresponding enhanced shell 
gap for the case of magic Z = 50. Fig.fTlTc) shows the 
behavior of the mid-shell region around A^ = 66 (exactly 
in between two shell closures: 50 and 82): the neutron 
shell gap for A^ = 66 is between the one for A^ = 65 and 
A^ = 67. Fig. fTlT b) presents the shell gap around A^ = 40. 
For A^ = 40 a strong difference (like for A^ = 82) is not 
visible and N = 40 is distinct from neither A^ — 39 nor 
41. Note that the A^ = 39 mid-shell gap is larger than 
those of A^ = 38 and 40 for several values of Z, especially 
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FIG. 9: Two-neutron separation energies (S2n) minus a lin- 
ear function of N around N — 'iG (left), and the strong shell 
closure N = 82 (right), for comparison. 



FIG. 10: (a) Shell gap as a function of the proton number Z 
for different magic and mid-shell neutron numbers N. N=16, 
28, 50, 82 correspond to shell closures, = 39 and 66 are 
exactly between two shell closures (called mid-shell), TV = 40 
is under investigation. Data are from fl^] . 



for Z = 28, unlike the N = 66 mid-shell behavior. This 
shows that N = 38, 39, and 40 do not have the behavior 
we could have expected from observation in other mass 
regions. However, in summary, no shell closure at TV = 40 
is observed. 



a peak that would indicate a shell closure, as shown in 
Fig.[l2tb) where shell closure at TV = 28, 50, and 82 are 
clearly visible. 



C. The pairing gap 

The pairing gap from the four-point formula [74 

A4(TV, Z) 



(-1) 



N 



B{N + 1) ~3B{N) 
+ 3B(^N -1) - B{N -2)](9) 



A4(TV,Z) 



was chosen to study the pairing-energy behavior. A peak 
is expected for magic numbers and a trough at mid-shell. 

The pairing gap as a function of neutron number is 
presented in Fig.^a) for Z = 28 - 32. At the TV = 
39 mid-shell, there is a trough for Z = 31 - but not for 
Z = 29. A similar behavior is seen at TV 66 (82-50 mid- 
shell). The odd-Z nucHdes have a lower pairing gap and 
while germanium {Z = 32) shows no particular structure, 
nickel (Z = 28) shows a strong mid-shell trough and not 



D. Comparison with the Bethe-Weizsacker formula 

The Bethe-Weizsacker formula was given in eq. (HJ . We 
adapt the version of Pearson [7^ , with a pairing term of 
Fletcher [tHI- Thus, the binding energy per nucleon is 
given by 



3e^ 



A 



5ro 



-4/3 



-\-apA~y- 



i/(-l)^ + (-l) 



N 



(10) 



with / = (TV — Z)/A. The parameters are ayoi = 
-15.65 MeV, a^/ = 17.63 MeV, a,, = -25.60 MeV 
which is the parameter of surface symmetry term intro- 



duced by Myers and Swiatecki [77 



sym 



27.72 MeV, 



ro = 1.233 fm with tq the constant used in the radius 
estimation R ~ tqA^/^, ap = —7 MeV the pairing term. 
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FIG. 11: Shell gap as a function of the proton number Z for a) 
= 81 — 83 with the N=82 magic number well distinguished 
from N=81 and 83, b) N = 2,9,- 41, and AT = 65 - 67 with 
N=66 representing a mid-shell number in between N=50 and 
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FIG. 12: (a) Pairing gap energy as a function of neutron 
number for the investigated elements: nickel, copper, and gal- 
lium, as well as zinc and germanium, (b) Pairing gap energy 
as a function of neutron number for Z=27-59. Shell closures 
at N=28, 50, and 82 are clearly visible, the N=66 mid-shell 
is indicated. 
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FIG. 13: Difference between the experimental mass values 
from this work and from AME2003 data and theoretical 
masses from the "Bethe- Weizsacker formula" as a function of 
proton number, for several magic neutron numbers and for 
N=40. 



FIG. 14: Difference between the masses predicted by the 
Bethe- Weizsacker formula (eq. \10]) and the experimental val- 
ues as a function of N for Z — 28, 29, and 31. Data are from 
this work complemented by [4^j- 



and y — 0.4. This formula contains no specific term for 
shell effects so the formula may not be a good way to pre- 
dict exotic mass values. However this makes it a "neutral" 
indicator for shell structures fsee [78|). 

To this end, the modified Weizsacker formula [eq. (fTO| ] 
is subtracted from known masses (divided by A) . The dif- 
ference between the experimental values and the formula 
clearly reveals the shell closures at TV = 28, 50, 82 and 
126, reaching up to 15 MeV for TV = 50 and TV = 82 (see 
Fig. 1 injli). 

Fig. [13] presents the difference between the experimen- 
tal results obtained from this work (complemented with 
AME2003 data) with the "Bethe- Weizsacker formula" 
[eq. fTO|) ] as a function of Z for various magic neutron 
numbers, including TV = 40. As with the shell gaps, the 
cases where N — Z show the strongest effects, as does the 
case of 5Q^Sn82. Interestingly enough, the case of 2iNi4o 
does show a dip of about 2 MeV, although only about 
20% the effect of ll^Sns2. 

When the difference in mass values is examined iso- 
topically as a function of neutron number fFig. fT4| , how- 
ever, there is no indication of a shell, or even sub-shell 
closure. The pseudo-parabolic behavior of the curve in 
Fig. HH shows some indentation around TV = 40 but noth- 



ing that we could claim to be "magic". 



VI. CONCLUSION 

The high-precision mass measurements performed at 
ISOLTRAP on over 30 short-lived neutron-rich isotopes 
of nickel, copper, and gallium have allowed us to rather 
finely study the mass surface - and its derivatives - 
around the interesting region of Z = 28 and TV = 40. 
No behavior resembling that of known magic numbers 
has been found, unlike the analog case oi Z = 40, where 
the TV = 56 sub-shell closure is visible. As much as an 
TV = 40 ((^5/2) sub-shell could exist, there is no clear in- 
dication for such a sub-shell closure from these measure- 
ments. While the pairing gap energy clearly indicates 
that there is no shell closure in this region, a competing 
mid-shell stabilization effect might be present. The com- 
parison with the Bethe- Weizsacker formula shows some 
fine structure around TV = 39,40 but no indication of 
the presence of a shell, or sub-shell closure. The shell 
gap evaluation shows anomalous behavior for TV = 39 as 
well as for TV = 40, perhaps due again to the competition 
between a sub-shell closure at 40 and the mid-shell at 39. 

Recalling again the words of Bohr and Mottelson, "it 
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is relatively difficult to discern the nuclear shell struc- 
ture as long as the main information on nuclei is con- 
fined to binding energies". While they are a necessary 
ingredient, it is not sufficient for explaining the prob- 
lem at hand since the bindin g en ergies are in opposi- 
tion with results on the -B(£'2') |la|. Thus, more detailed 
spectroscopy measurements, including the 5— factor, as 
suggested by Langanke et al. [l^ , and more theoretical 
work, are called for to understand the various phenomena 
arising from mass-surface studies. 
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